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ABSTRACT 

Aims. We present dense grids of stellar models suitable for comparison with observable quantities measured with great precision, 
such as those derived from binary systems or planet-hosting stars. 

Methods. We computed new Geneva models without rotation at metallicities Z = 0.006, 0.01, 0.014, 0.02, 0.03 and 0.04 (i.e. [Fe/H] 
from -0.33 to +0.54) and with mass in small steps from 0.5 to 3.5 M G . Great care was taken in the procedure for interpolating between 
tracks in order to compute isochrones. 

Results. Several properties of our grids are presented as a function of stellar mass and metallicity. Those include surface properties 
in the Hertzsprung-Russell diagram, internal properties including mean stellar density, sizes of the convective cores, and global 
asteroseismic properties. 

Conclusions. We checked our interpolation procedure and compared interpolated tracks with computed tracks. The deviations are less 
than 1 % in radius and effective temperatures for most of the cases considered. We also checked that the present isochrones provide nice 
fits to four couples of observed detached binaries and to the observed sequences of the open clusters NGC 3532 and M67. Including 
atomic diffusion in our models with M< 1.1 M Q leads to variations in the surface abundances that should be taken into account when 
comparing with observational data of stars with measured metallicities. For that purpose, iso-Z surf lines are computed. These can be 
requested for download from a dedicated web page together with tracks at masses and metallicities within the limits covered by the 
grids. The validity of the relations linking Z and [Fe/H] is also re-assessed in light of the surface abundance variations in low-mass 
stars. 
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1. Introduction 

Accurate grids of stellar models at various metallicities remain 
a basic tool in many areas of astrophysics, such as the study of 
stellar clusters, of stellar populations, of stellar nucleosynthesis, 
or of chemical evolution in galaxies. Two observational develop- 
ments have further strengthened the need for ever more accurate 
stellar model predictions in recent years. The first concerns the 
improvements in stellar parameters determination from observa- 
tions. We can mention in this respect the determination to better 
than 3% of stellar masses and radii of almost two hundred stars 



in eclipsing binary systems (Torres et al. 2010 1 or the precise 
measurement of the characteristics of planet-hosting stars. The 
second development is related to the advent of large-scale sur- 
veys, initiated in the nineties with the EROS, MACHO, OGLE, 
Hipparcos, or SDSS experiments, to cite only a few. Those sur- 
veys provide statistically significant data on millions to tens of 
millions of stars, and future ones are expected to increase the 
number of stars by one or two orders of magnitude. We cite 
in this respect the future Gaia mission, planned for a launch 
in 2013, which will monitor 1 billion stars in our Galaxy, with 
distance determinations to better than 1 % for about 20 million 
of them allowing the building of the Hertzsprung-Russell (HR) 
diagram to an unprecedented degree of precision. The Large 
Scale Synoptic Survey is another example of a future large scale- 
survey, which will complement Gaia by, among other goals, ob- 



serving fainter objects. The start of its operation is planned for 
the end of this decade. 

Significant improvements have been made to the field of stel- 
lar modeling since the first extensive grids of Geneva stellar 



models published about two decades ago ( S dialler et al.|l992 
and subsequent papers). They concern the equation of state, 
atomic diffusion, rotation or the solar abundances, to cite only 
a few. The impact of their inclusion in our stellar models has 
been presented in a few papers (e.g. Eggenberger et al. 2010a|b| l, 
justifying the need to update our previous grids of stellar models. 

New grids of Geneva models covering stellar masses from 
0.8 to 120 M at solar metallicity and computed with and with- 
out rotation are presented in a separate paper (Ekstrom et al. 
|2011[ ). In the present paper, we present dense grids of mod- 
els covering stellar masses from 0.5 to 3.5 M and metallicities 
from Z = 0.006 to 0.04 with a distribution of points in the mass- 
metallicity plane dense enough to allow an accurate interpolation 
between computed models, but restricted to models without ro- 
tation and to phases from the zero-age main sequence (ZAMS) 
up to the base of the red giant branch (RGB). 

The originality of the present study resides mainly in three 
aspects. First, we use the solar abundances given by Asplund et 
al. (2005) as initial abundances. Second, we provide tracks for a 
very dense sampling in initial masses for six metallicities. Third, 
we construct an interpolation procedure between the tracks, al- 
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Table 1. Initial chemical abundances of our grids. 
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[Fe/H] wt 


■^init 




•Zinit/^init 


0.006 


-0.331 


0.7383 


0.2557 


0.0081 


0.010 


-0.103 


0.7291 


0.2609 


0.0137 


0.014 


+0.048 


0.7200 


0.2660 


0.0194 


0.020 


+0.212 


0.7063 


0.2737 


0.0283 


0.030 


+0.402 


0.6834 


0.2866 


0.0439 


0.040 


+0.542 


0.6606 


0.2994 


0.0606 



manner than models using older solar abundances do (see e.g. 
|Basu & An tia 2008 and references therein). 

Atomic diffusion Atomic diffusion due to concentration and 
thermal gradients is included in the computation of models with 
initial masses below 1.1 M Q , following the Chapman-Enskog 



lowing us to efficiently handle the change in the morphology of 
the tracks when the mass varies. 

The paper is divided as follows. Secfion[2]presents the ingre- 
dients of the stellar evolution models and Sect. [3] some general 
properties of the grids. Section |4]describes the different types of 
tracks and isochrones constructed from the raw models. A com- 
parison of our model predictions with some observations and 
with some other tracks available in the literature is presented in 
Sect. |5] The data made available for download are described in 
Sect.|6] Conclusions are then drawn in Sect. [7] 

Three appendices complete the main body of the paper. 
Appendix [A] discusses the relation between the two standard 
representations of the metallicity, i.e. Z, the mass fraction of 
all elements heavier than helium, and [Fe/H], the logarithm of 
the iron to hydrogen abundance ratio relative to the solar ra- 
tio. Appendix [B] describes the procedure used to construct basic 
tracks, which consists in a well defined series of stellar mod- 
els chosen to insure easy interpolations in mass and metallicity. 
Appendix [C] describes the construction of isochrones. 



2. Ingredients of the stellar models 

The input physics included in the Geneva evolution code is de- 
scribed in detail in Eggenberger et al. (2008 ), while the chem- 
istry and nucleosynthesis used for the grids is described in 
Ekstrom et al. ( 201 1) . We highlight here only a few points rele- 



vant to the present dense grids. 



Chemical abundances The adopted initial chemical abun- 
dances are indicated in Table [T] (see Appendix [A] for more de- 
tails). For the mixture of the heavy elements, we adopt the ones 

( 2005} , except for Ne which we take from 3. Grids properties 



method ( Chapman & Cowling|1970 i. We refer to Eggenberger 
|et al.] (2008) for more details about the modeling of atomic dif- 
fusion in the Geneva evolution code. For masses above about 
1 . 1 M , helium and heavy elements are rapidly drained out of the 
envelope when only atomic diffusion due to concentration and 
thermal gradients is included to model those stars (which have 
a thin convective envelope). This results in surface abundances 
that are too low compared to observations. Indeed, in those stars, 
the effects of atomic diffusion due to radiative forces and/or to 
macroscopic turbulent transport become significant. These ef- 
fects mainly counteract those of atomic diffusion due to concen- 
tration and thermal gradients. Therefore, we assume here that 
above the (somewhat arbitrary) mass limit of 1.1 M , those mod- 
els with no atomic diffusion probably better fit the observations. 
We thus do not include atomic diffusion in stars with a mass 
equal to or greater than 1.1 M . 

Equation of state Models with M > 1.1 M are computed with 
the general equation of state (EOS) used in the Geneva evolution 
code that is well suited to massive stars (Schaller et al.|1992) . 
The computation of reliable models of solar-type stars requires, 
however, a more specific equation of state, because non ideal 
effects such as Coulomb interactions become important for low- 
mass stars. The OPAL EOS ( |Rogers & Nayfonov|2002l l is then 
used for models with initial masses below 1.1 M . We checked 
that the transition between models with diffusion and the OPAL 
EOS (M < 1 . 1 M ) and models without diffusion and with the 
general EOS (M> 1.1 M ) is smooth in the physical quantities 
at the center of our models. 



Mass loss No mass loss is applied to our models, which all 
have M<3.5M and are computed only up to the base of the 
RGB. 



from Asplund et al 
|Cunhaetal.|p006^ 



Convection and overshooting Convective zones are deter- 
mined with the Schwarzschild criterion. The convective core is 
extended with an overshoot parameter a m = d m /Hp of 0.10 for 
M > 1.7 M . Half of this overshooting, a ov = 0.05 is applied 
between 1.25 and 1.7 M , where the extent of the convective 
core, when it exists, is small. If d ow exceeds the dimension of 
the Schwarzschild core /?sch, then the total extension of the con- 
vective core is taken equal to/?schX(l+fl?ovei7#p)- This procedure 
avoids having an extension of the core superior to the radius of 
the initial core. No overshooting is applied below 1.25 M . 

The outer convective zone is treated according to the mixing 
length theory (Bohm-Vitense 1958 ), with a solar calibrated value 
of the mixing-length parameter otmlt = I /Hp equal to 1.6467 
based on models including atomic diffusion (see below). This 
low value of the mixing-length parameter directly results from 
the use of the new solar abundances from Asplund et al. (2005). 
It is, however, known that these abundances lead to solar mod- 
els that reproduce helioseismic constraints in a less satisfactory 



Our dense grids of stellar models are computed for initial metal- 
licities Z; n i t and masses M equal to 



0.006, 0.010, 0.014, 0.020, 0.030, 0.040 



M 



= 0.5 [0.1] 1.1 
1.1 [0.02] 1.3 
1.3 [0.1] 3.5 



M , 
M , 



(1) 



where the notation M a [dM] Mb stands for a range of masses be- 
tween M a and M\, in steps of dM. Table[T]gives some information 
on the initial abundances at each metallicity. 

Section |3.1| summarizes the model predictions in two dia- 
grams relating quantities derivable from observations, namely 
the classical HR diagram and the logT^ - log p/p Q diagram 
plotting the mean stellar density against effective temperature. 
Sectionl3~l2"1fhen addresses the onset of convection in the core of 



intermediate-mass stars, Sect. 3.3 discusses the ages of our mod- 
els, and Sect. 3.4 briefly presents the effect of atomic diffusion 
in low-mass models. 
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Fig. 1. Evolutionary tracks in the HR diagram at Z; n j t = 0.014. 
The main sequence is shown in thicker lines. Reference points 
defined in Sect. |4.1| are indicated by diamonds (ZAMS and 
TAMS) and filled circles (three reference points within the MS 
and one at the bottom of the RGB). For clarity, the tracks at stel- 
lar masses between 1.1 and 1.2 M Q and between 1.2 and 1.3 M Q 
are not shown. They are displayed in Fig.|4]for Zj n ; t = 0.006 and 
0.04. 
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log Teff 
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Fig. 2. Evolutionary tracks at Z; n i t = 0.006 (in blue, hot tracks) 
and 0.040 (in red, cool tracks). Only selected stellar masses are 
plotted for clarity, as labeled next to the ZAMS location of the 
low-metallicity grid. The main sequences are drawn in thin con- 
tinuous lines and the post-main sequences in thin dotted lines. 
Thick continuous lines denote the location of the ZAMS (hot 
sides) and of the onset of the RGB (cool sides). Also shown for 
each of the two metallicity grids are the isochrones at 13.5 Gyr 
(thick dashed lines on the low-luminosity side). 



3.1. Summary diagrams 

Hertzsprung-Russell diagram The HR diagram of our 
Zinit = 0.014 tracks is shown in Fig.[T] with the ZAMS, the termi- 
nal age main sequence (TAMS), and the base of the RGB. The 
base of the RGB is defined as the point on the post-MS track 
where the slope of the line connecting the TAMS and the running 
point on the post-MS tracks starts to increase (see Sect. B.2i. 
Three other reference points chosen on the MS are defined in 
Sect. |4. l| and used for constructing the basic tracks. Two features 
are worth noticing in the HR diagram. First, a hook appears at the 
end of the MS phase in the morphology of the MS tracks at stel- 
lar masses above 1.1 M Q . The filled red and green points locate 
respectively the red and the subsequent blue turn off points. The 
hooks lead to a degeneracy in the HR diagram, whereby a given 
set (L, Trff) at a given metallicity can correspond to up to three 
stellar models, each with a different stellar mass. Second, from 
ZAMS to the red turn-off point (or to the TAMS in the absence 
of hook), 7^.ff strictly increases with time for masses below 1 M Q , 
strictly decreases with time for stellar masses above 1.4 M G , and 
has a decreasing-then-increasing behavior for masses in between 
those limits. 



The area covered by all our models in the HR diagram is 
shown in Fig. [2] from Zj n ; t = 0.006 (hotter tracks) to 0.04 (cooler 
tracks). The isochrones at 13.5 Gyr plotted in the figure start 
on the 0.5 M Q tracks and remain very close to the ZAMS up 
to a mass around 0.75 M because very low-mass stars do not 
evolve much in the time span of 13 Gyr (the MS lifetime of the 
0.5 M Q star ranges between 69 Gyr at Zj„; t = 0.006 and 76 Gyr at 
Zinit = 0.04, much greater than the age of the Universe). 

The metallicity dependence of the tracks is shown in Fig. [3] 
for three stellar masses. The morphology of the tracks at 1.2 M Q 
is particularly interesting as it shows a metallicity dependence 
of the MS hook in the HR diagram due to the sensitivity of the 
onset of convection to metallicity (see Sect. 3.2 1. The tracks in 
the critical stellar mass range between 1.1 and 1.3 M Q , where 
convection sets in in the core of MS stars, are shown in more 
detail in Fig. |4]for the two extreme metallicities considered in 
our grids (see discussion in Sect. 3.2 1. 



Effective temperature - mean density diagram The mean den- 
sity p = M/( ^nR 3 ) of a star can be determined with relatively 
good accuracy in systems hosting exoplanets (e.g. |Sozzetti et al.] 
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Fig. 3. Evolutionary tracks of the 0.7, 1.2, and 2M models at 
the different metallicities of our grids as labeled in the figure. 
Reference points are indicated on the tracks by filled diamonds 
and circles in a similar way to Fig. [T] 
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Fig. 4. Evolutionary tracks at Z; n j t = 0.04 (upper panel) and 0.006 
(lower panel) for the intermediate-mass stars from 1.1 to 1.3 M Q , 
as labeled in the panels. Reference points are indicated on the 
tracks by filled diamonds and circles in a similar way as in Fig.[T] 



2007 Gillon et al.|20 09). This quantity generally decreases with 
age for a given star, as shown in Fig. [5] which displays the tracks 
followed by the Zi n ; t = 0.014 models in the logoff - log p/p G 
plane. 

The range of stellar mean densities and effective tempera- 
tures covered by our grids is shown in Fig. [6] for the two ex- 
treme metallicities considered here. For the very low-mass stars 
(M<0.8), however, as in the HR diagram, a portion of the 
\ogT e ff - log p/p Q plane is unreachable within the age of the 
Universe. The frontier of the realistic region is indicated in Fig. [6] 
by the 13.5 Gyr isochrones, starting on the 0.5 M Q tracks close 
to the ZAMS and ending at the base of the RGB on the 0.9 (at 



Z init = 0.006) and ~1 M Q (at Z init = 0.04) tracks. 



3.2. Onset of convection in the core 

There is a mass limit below which stellar models burn hydro- 
gen radiatively throughout the MS phase, and another mass limit 
above which the core is convective all through the MS phase. 
Those mass limits depend on metallicity. For Zinit = 0.014 stars, 
they are ~1.20M Q and ~1.24M , respectively. Between those 
two mass limits, a star may have a convective core at ZAMS 
that disappears with age and/or a convective core that appears 
in the middle of the MS, grows thereafter and disappears at the 
end of the MS. This is seen in Fig. [7] which displays the evo- 
lution of the mass of the convective core in the 1.1 to 1 .26 M Q 
models as a function of relative time on the MS at three dif- 
ferent metallicities. In the 1.20 M Q model at Z^, = 0.014 (mid- 
dle panel of the figure), for example, a convective core exists at 
ZAMS and disappears after 0.74 Gyr (when the hydrogen mass 



fraction at the center of the star, X c , equals 0.61). A convective 
core then reappears at 1.89 Gyr (X c =0.22), reaches its maxi- 
mum growth in mass at 3.79 Gyr (X c =0.13), and disappears at 
4.18 Gyr (X c = 0.015). The MS ends at 4.93 Gyr (X c = 10~ 5 ). In 
the 1.22 M Q stellar model (MS duration of 4.62 Gyr), the initial 
convective core disappears at an age of 1.15 Gyr (X c = 0.53), and 
reappears at 2.25 Gyr (X c = 0.32). 

There are two processes at work in making the core unstable 
against convection. The first process is due to the fusion of 3 He 
through the reaction 3 He + 3 He — > 4 He + 2 p. It is responsible for 
the occurrence of a convective core right from the ZAMS in stars 
with M > 1 . 1 M G (at Z illit = 0.014). That the 3 He + 3 He reaction is 
very energetic (it depends on the 18 th power of temperature at 
12 x 10 6 K) leads in those stars to the appearance of a convec- 
tive core at ZAMS, with a mass that decreases with time as 3 He 
burns. The core temperatures at ZAMS are shown in Fig. [8] The 
lower stellar mass at which the core is convective on the ZAMS 
corresponds to the beginning of the thick line at the left. It is 
computed by extrapolation from the mass of the convective core 
as a function of stellar mass, and is equal to 0.984, 1.055, 1.113, 
1.176, 1.215 and 1.236 M B at Z init = 0.006, 0.010, 0.014, 0.020, 
0.030, and 0.040, respectively. 

Thus one sees that the minimum mass for having a con- 
vective core on the ZAMS increases when the metallicity in- 
creases. This comes from the fact that, for a given mass, hydro- 
gen burning occurs at higher temperature when the metallicity 
decreases (since the luminosity of the star is higher). This will 
allow the 3 He + 3 He reaction to release a greater amount of en- 
ergy since the initial abundance of 3 He does not depend (or in 
a very marginal way, see e.g. |Lagarde et al.|2Ql 1[ > on the initial 



4 



N. Mowlavi et al.: Stellar mass and age determinations 




3,9 3,8 
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Fig. 5. Evolution in the log r e ff-logp/p Q diagram at Z; n j t = 0.014. 
See Fig.[T]for the explanation of the symbols. 
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Fig. 6. Evolution in the log 7^ - log p/p e diagram at Zinit = 0.006 
and 0.04. See Fig.[2]for the explanation of the symbols. 



metallicity. This favors the appearance of a convective core at 
lower initial mass at low metallicity. 

The second process that makes the central regions of a main 
sequence star unstable against convection is the high tempera- 
ture dependence of the CNO burning. It comes into play when 
the temperature reaches a level high enough for the CNO cy- 
cles to outweigh the p-p chains energetically. This occurs for 
stellar masses higher than about 1 .2 M Q , and is at the origin of 
the growth of the convective core in the second part of the MS 
as shown in Fig. [7] The evolution of the convective core in the 
models up to 2 M Q is shown in Fig. [9] It is only for stellar masses 
above ~1.7 M that the mass of the convective core is a strictly 
decreasing function of time during the whole MS, since the CNO 
chains are highly efficient from the ZAMS in those stars. 

The mass of the convective core is therefore not always a 
strictly decreasing function of time during the MS. The growth 
of the convective core with time, especially at stellar masses be- 
tween 1.2 and 1.25 M Q , implies mixing of fresh hydrogen from 
the outer layers into the H-depleting core. The discrete time 
steps needed to compute evolutionary models result in sudden 
increases in the hydrogen abundance in the core, leading to dis- 
continuous abundance profiles at its border. These sudden ad- 
ditions of a discrete amount of hydrogen to the core may lead 



We add two remarks concerning convective cores in MS 
models. The first remark concerns the sensitivity of their extent 



at ZAMS to pre-MS modeling. [Morel et al. (2000 1 performed 



a study of solar models and compared the evolution of models 
when the pre-MS phase is accounted for and when it is not. The 
authors conclude that both sets of models become nearly identi- 
cal as soon as the nuclear reactions start to work at equilibrium. 
With the ZAMS defined when Z C (H) has decreased by 0.25% 
in our models (see Sect. |4. 1[ ), a value chosen to ensure that nu- 
clear equilibrium is reached for all masses and metallicities in 
our grids, our MS models are independent of the pre-MS evo- 
lution. We also note that with our definition of the ZAMS, our 
solar track does not have a convective core at ZAMS. The con- 
vective core that develops at the early stage of the solar track 
(see |Morel et al.|2000l ) has exhausted by the time of the ZAMS. 

The second remark concerns the sensitivity of the extent of 
the convective cores on the constitutive physics of the mod- 
els. [Magic~et~ar (2010) show, for example, how the constitu- 



tive physics of the models, such as the extension of the con- 
vective core, can reproduce a hook in the critical mass range 
1.2 — 1.3 M Q . The study was done in relation to the open clus- 
ter M67, which has a turn-off mass precisely in that mass range 
and which displays a clear hook in the color-magnitude diagram 
(CMD). We discuss this case further in Sect. [5] 



to 'spikes' in the hydrogen mass fraction with time. Figure 10 

which displays the evolution of X c with time, shows that such ^ <j Main sequence lifetimes 



spikes occur for stellar models in the range 1.20 - 1.24M . The 
hydrogen mass fraction in the core of MS stars is thus not nec- 
essarily a strictly decreasing function of time either. 



Figure[TT]displays the MS lifetimes (upper panel), fn, as a func- 
tion of stellar mass at three metallicities. It is a well known de- 
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Fig. 7. Evolution of the mass of the convective core as a func- 
tion of the relative time spent on the main sequence (0 is ZAMS 
and 1 is TAMS) for models at Z init = 0.040 (upper panel), 0.014 
(middle panel), and 0.006 (lower panel). The models are for stel- 
lar masses of, from top to bottom, 1.26, 1.24, 1.22, and 1.20 M Q 
(as labeled on the continuous lines in all panels), 1.18 M (long 
dashed lines in the middle and lower panels), and 1.1OM (short 
dashed line in the lower panel). 



creasing function of stellar mass, mainly due to the increase in 
stellar luminosity with stellar mass. The time needed to reach the 
base of the RGB from the TAMS is plotted in the lower panel of 
Fig. [TT] It amounts to about 30% of ?h for models with masses 
below ~1.2M , decreases rapidly above that mass to ~3% at 
1.8 M , and then does so more slowly for more massive stars, 
reaching -2% at 3.5 M Q . 

The sensitivity of f H with stellar mass can be analyzed, 
at any given metallicity, by comparison with the lifetimes at 
Zinit = 0.014. This comparison, shown in Fig. 12 reveals a max- 



imum sensitivity of ?h to Z at stellar masses between 0.8 and 
1.5 M . At lower masses, tu becomes less sensitive to Z, since 
it is almost equal in all 0.5 M models with Zj mt > Z 0> j mt . This is 
due to the independence of the p-p chain to metals and to the 
lower sensitivity of the p-p reaction rates to temperature. 

The sensitivity of ?h with Z is shown in Fig . [T3| The plot 
confirms the results highlighted above from Fig. |12| according 
to which the highest sensitivity is around 1 .2 M Q and the low- 
est at the smallest mass 0.5 M Q considered in our grids. The 
figure reveals, however, a pattern of ?h as a function of Z that 
differs from the pattern found in Mowlavi et al. ( 1998 1 (here- 



after MMM98) with the old Geneva models. While the 3 M 
was found by MMM98 to indeed be less sensitive to metallic- 
ity than the 1 M Q models (see their Fig. 5), in agreement with 
the results of this paper, the models at metallicities higher than 
solar showed a trend of decreasing ?h with Z, which is not found 
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Fig. 8. Upper panel: Central temperature, in million Kelvin, at 
ZAMS of stellar models at different metallicities as a function 
of the stellar mass of the models. Each dot represents a com- 
puted model. The thick lines indicate the range of stellar masses 
for each metallicity, where a convective core is present on the 
ZAMS. The dotted line connects the core temperature at the low- 
est stellar mass having a convective core, for the different metal- 
licities. Lower panel: Same as upper panel, but for the central 
density. 
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Fig. 9. Same as Fig. [7] but for stellar masses of, from bottom to 
top, 1.3 to 2 M in steps of 0.1 M as labeled on the lines. 
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Fig. 10. Same as Fig. [7] but for the evolution of the hydrogen 
mass fraction at the center of the star in the models from 1.20 to 
1.26 M B . 



in this paper. This results first from the adoption in the present 
work of the solar metallicity of Asplund et al. (2005), implying 
that a given value of Zj n i t /Z 0; j n i t corresponds now to a lower value 
of Zi n i t . In the present work, we therefore never reach the very 
extreme metallicities explored in MMM98. Second, the slope 
AY/AZ- 1.2857 adopted in this paper is a factor of almost two 
lower than in MMM98. This prevents to have the significant 
drop of the initial H abundance recorded at high metallicities 
in MMM98 and which contributed to their decrease of ?h with 
metallicity. Moreover, many physical ingredients have been up- 
dated, such as the opacities, some nuclear reaction rates, and the 
overshooting parameter, which may of course contribute to ex- 
plaining this difference as well. 

In the cases where the mass of the convective core grows 
with time, the addition of fresh hydrogen naturally increases the 
main sequence lifetime. The question then arises whether nu- 
merical simulations do adequately model this increase. The evo- 
lution of the mass of the convective core with time (Fig. |7]i sug- 
gests that the global evolution of the convective core during the 
MS is not affected much by the numerical spikes. It must also be 
noted that the effective temperature and the surface luminosity 
of the models are evolving smoothly with time irrespective of 
the occurrence of spikes. We thus assume that the MS lifetimes 
are not affected by those spikes. The situation is quite different 
from what happens during the core helium-burning phase where 
spikes in the core helium abundance appear in the growing core 
at a late stage of helium burning, when the mass fraction of he- 
lium has decreased below 10%. The spikes at that stage lead to a 
significant increase in the core helium-burning lifetime, the later 
becoming sensitive to the numerical parameters adopted in the 
simulations. We do not expect that difficulty during the MS. 
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Fig. 11. Upper panel: Main sequence lifetimes in Gyr as a 
function of initial stellar mass at three different metallicities: 
Zinit = 0.006 (short-dashed blue line), 0.014 (continuous black 
line), and 0.040 (long-dashed red line). Lower panel: Same as 
upper panel, but for the time, relative to the MS lifetime, needed 
to the models to cross the HR diagram from the end of the MS 
to the bottom of the RGB. 
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Fig. 12. Main sequence lifetimes, as a function of initial stellar 
mass, of models at different metallicities as labeled next to the 
curves. The lifetimes are displayed relative to the lifetime of the 
considered stellar mass at Zi n ; t = 0.014. 



3.4. Surface abundances 

Atomic diffusion leads to nonconstant surface abundances with 
time. In the 1.0 M Q , Zi n ; t = 0.014 model, for example, the sur- 
face hydrogen mass fraction X sur f increases by about 7% during 
the main sequence (MS), while that of helium, T sur f , decreases 
by almost 20%. This is displayed in, respectively, the upper and 
second panels from the top of Fig. 14 As a result, the surface 
metallicity decreases by 1 1 % during the MS of that star, reach- 
ing Z sur f = 0.0124 at the end of that phase (third panel from top 
of Fig. 



14) 



It evolves back to its initial value during the post- 
MS phase as the convective envelope deepens into the star and 
homogenizes back the diffused elements. Similar behaviors are 
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Fig. 13. Upper panel: Main sequence lifetimes as a function of 
the logarithm of the metallicity relative to solar, for models from 
0.5 to 1 . 1 M (bottom to top curves) as labeled on the right of 
the curves. The lifetimes are normalized for each stellar mass to 
its value at Z; n i t = 0.014. Lower panel: Same as upper panel, but 
for selected models from 1.1 to 3.5 M G (top to bottom curves). 



seen in the evolution of the surface abundances at stellar masses 
below 1.0 M Q , 



14 for the 0.9 and 0.8 M Pj 



as confirmed in Fig. 
models at Z; n ; t = 0.014. The change of the surface metallicity due 
to diffusion is thus not negligible during the MS phase. 



Figure 15 further displays the evolution of the surface metal- 
licity during the MS for all our models that include diffusion at 
Zinit = 0.04 (upper panel), 0.014 (middle panel) and 0.006 (low- 
est panel). The fraction of metal that depletes from the surface 
owing to atomic diffusion into the deeper layers of the star de- 
creases with increasing metallicity for a given stellar mass, due 
to the different densities at the base of the convective envelope. 
As a result, 10 % of the metal mass fraction is depleted by the 



end of the MS in the 
0.006. 



1 M model at Z init = 0.04, and 20 % at 



3.5. Global asteroseismic quantities 

By providing a wealth of information on the internal struc- 
ture of the Sun, the solar five-minute oscillations stimulated 
various attempts to obtain similar observations for other stars. 
These past years, a growing number of detections of solar-like 
oscillations have been obtained for different kinds of stars 
by ground-based spectrographs and space missions. Scaling 
relations relating asteroseismic parameters to global stellar 
properties are largely used to analyze these data. These global 
asteroseismic quantitites are indeed of great help to constrain 
stellar parameters and to obtain thereby new information about 
stellar evolution without the need to perform a full asteroseismic 
analysis (see e.g. |Stello et al.|200"9"b"}|Migho et al.|2009t|Chaplin| 
et al. 201 lb l. In this section, we provide some asteroseismic 
properties of our models that can be directly computed from 
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Fig. 14. Hydrogen mass fraction (upper panel), helium mass 
fraction (second panel from top) and metallicity (two lower pan- 
els) at the surface of the 0.8 (dashed-short), 0.9 (dashed-long) 
and 1.0 M G (continuous) Zj n ; t = 0.014 models as a function of 
age. Filled diamonds indicate the end of the main sequence 
phase. Filled circles on the 1 M Q tracks locate the solar model 
at an age of 4.57 Gyr. 



their global properties using scaling relations: 



the large frequency separation Av that is proportional to the 
square root of the mean density of the star (e.g. Ulrich 1986 ): 



Av = Av x 



(M/M )' 



0.5 



(2) 



(R/R e y- 5 ' 

with Av = 134.9 fjHz ( |Chaplin et al.|2011a) ; 

- the frequency v max at the peak of the power envelope of the 
oscillations where the modes exhibit their strongest ampli- 
tudes, given by (Brown et al. 1991 Kjeldsen & Bedding 
[1995) : 



(«/«o) Ueff/4ff,o) 



(3) 



with v max , = 3150 //Hz (Cha plin et al.|201 la[ ); 



- and the maximum oscillation amplitude A max (e.g. Chaplin 
|et al.|2011a) : 



M 



?eff 



'eff.o 



(4) 
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Fig. 15. Evolution of the surface metallicity during the MS in 
all our models with M < 1 M Q at Z init = 0.04 (upper panel), 0.014 
(middle panel), and 0.006 (lowest panel). The X-axis has been 
restricted to the age interval [0-20] Gyr. 




I 0,0000 
S> -0,0025 
-0,0050 



2,0 2,5 3,0 3,5 4,0 
log(nu_max/microHz) 



45 



Fig. 16. Top panel: Evolutionary tracks at Z m n- 0.014 in the 
Av/v max versus log(v max ) diagram. For clarity, only the 0.5, 1.1, 
and 3.5 M Q tracks are plotted, as labeled next to the tracks. The 
main sequence is shown in thicker lines. The ZAMS is located 
on the right of each track. Reference points are indicated with 
symbols as in Fig. [TJ The dashed line is a straight line passing 
through the ZAMS and the first MS hook of the 3.5 M Q star, 
given by Eq. [5] The frequencies are expressed in /iHz. Bottom 
panel: same as top panel, but for the residual of Av/v max with 
respect to the straight line in the top panel. 
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Fig. 17. Same as Fig. 



16 



but for log(Av) versus log(v max ). The 
dashed line in the top panel is a straight line with slope 0.77 
passing through the ZAMS of the 1.1 M Q model. 
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The 



dashed line in the top panel is a straight line passing through the 
ZAMS and TAMS of the 0.5 M G star, given by Eq. [6] 



where the solar value A max0 of the maximum oscilla- 
tion amplitude depends on the instrument used to mea- 
sure the oscillations. For the Kepler mission, for example 
A 



Mnax.O 



2.5 ppm (Chaplin et al. 201 la I. Following Kjeldsen 
2 and take 5 = 0.838 and 



& Bedding ( |1995|l we adopt r - 
t = 1.32 from |Huber etaLl l T^rT) . 
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We do not test here the validity of the scaling relations |2]to|4] but 
simply assume that these relations hold for the different masses 
and metallicities computed for these grids. 

The asteroseismic data of our grids are summarized in the 
top panels of Figs. 16 to 18 Figure [redisplays Av/ 
log Av and Fig 



18 



, Fig. 17 



jAmaxMmax.o, all three versus log(v max ). The 
oscillation frequencies and the large frequency separation de- 
crease during the evolution of the stars, as expected from the de- 
crease inof their mean density (cf. Fig. [5]), while the amplitude 
of their oscillation increases. 

The MS tracks are seen to obey an almost linear relation in 
each of the three figures. Suggested linear relations are shown 
by the dashed lines in the top panels of the figures. The residuals 
of the tracks with respect to that linear relation are shown in the 
bottom panels. In the [log(v max ), Av/v max ] plane (Fig. [IS), the 
linear relation is constructed based on the ZAMS and the point 
of first hook of the 3.5 M Q star. It is given by 



Av / v max = 0.128 - 0.0241 log(v max ) . 



(5) 



The coefficients of this relation vary only slightly with metal- 
licity, the relation being Av / v max = 0.133 - 0.0249 log(v max ) 
at Zjnit = 0.006 and Av/v max = 0.124 - 0.0235 log(v max ) at 
Zinit =0.04. 

In the [log(v max ), log(Av)] plane (Fig. [XT}, we adopt a linear 
relation with a slope of 0.77 as derived observationally (Stello 
et al.|2009a) and passing through the ZAMS of the 1.1 M Q track. 
The residuals displayed in the bottom panel of Fig. [l7]show very 
good agreement with our predictions. Our models actually pre- 
dict a slope of 0.74 for the MS track of the 0.5 and 1.1 M Q , and 
a larger slope of 0.82 for the MS track of the 3.5 M . 

In the [log(v max ), A max /A maxG ] plane, the linear relation is 
constructed based on the ZAMS and the TAMS of the 0.5 M Q . It 
is given by 



log(A max /A max , ) = 3.072 - 0.8734 log(v max ) . 



(6) 



The coefficients of this relation vary only slightly with metal- 
licity as well. The relation is log(A max /A max0 ) = 3.111 - 
0.8789 log(v max ) at Z init = 0.006 and log(A max /A max , ) = 3.037- 
0.8704 logOw) at Z m = 0.04. 

4. From tracks to isochrones 

We computed 234 raw tracks of evolutionary model, one for 
each {Zinit, M) couple taken from the list Q. The raw tracks 
start with an initial model close to the ZAMS and stop when 
the stellar model reaches the base of the RGB. The time steps 
between two successive models in each sequence are chosen to 
ensure good modeling of both the chemical (nucleosynthesis and 
atomic diffusion) and physical evolution. They vary with mass, 
metallicity and stage of evolution of the stellar models. 

In practice, for computing isochrones, we use tracks, here- 
after called basic tracks, described by fewer models than those 
defining the raw tracks. The models in the basic tracks are cho- 
sen in order to ensure reliable interpolation in mass or in metal- 
licity between them. Interpolations are necessitated because, in 
order to compute isochrones, we need to construct tracks for ini- 
tial masses (and metallicities) that are not among the 234 raw 
tracks. Besides this, these interpolations have to be made be- 
tween models that are at equivalent evolutionary stages. By this, 
we mean models that are either at the same evolutionary stage 
(for instance, the same mass fraction of hydrogen in the center 
or at the same relative time) or at the same position in the HR 
diagram (for instance, at the red turn off). We construct in this 



way 234 basic tracks, each consisting of 601 well chosen stellar 
models. In each basic track, the models with the same numbering 
are all equivalent, in the sense defined above. Each basic track is 
identified by a combination {Zj n ;t, M) taken from ([T]i. 

As explained above, evolutionary tracks at a {Z; n it, M] couple 
not available in the set of basic tracks can be computed by inter- 
polating among basic tracks. They are thus not stored in files, but 
can be computed on demand through the web page (see below). 
They are called interpolated tracks. 

In addition to the evolutionary tracks, we also construct lines 
in the HR diagram linking models of equal initial mass, having 
different initial metallicities and ages, but showing equal surface 
abundances. We call these lines iso-Z sur f lines. Each iso-Z sur f line 
is identified by a pair {Z sur f, M). 

Finally, to be complete, we recall that isochrones are lines 
linking models having the same age and initial metallicity. They 
are identified by a couple {Z; n it, age). The basic tracks, interpo- 
lated tracks, iso-Z sul f lines and isochrones are respectively de- 
scribed in Sects. |4~Tl |4~2l |431 and|474l 



4.1. Basic tracks 

Basic evolutionary tracks are constructed from the raw evolu- 
tionary tracks by first identifying reference points on each track 
in the log R - log L diagram, and then defining 'equivalent' mod- 
els by regularly distributing them between the reference points. 
We refer to Appendix [B] for a detailed description of the refer- 
ence points. 

The ZAMS is defined at the time when X C (H) has decreased 
by 0.25% at the center of the star. It ensures that the models have 
reached a gravitationally stable configuration. This is shown in 
Fig. [19] which displays, among other quantities, the evolution of 
the stellar radius (second panel from bottom) and of the central 
temperature (third panel from bottom) as a function of the frac- 
tion of central H burned in the case of Z^ = 0.04 models. The 
equilibrium state is reached at an even earlier stage of H burning 
in the grids of lower metallicities. 

The time taken by our initial models to reach the ZAMS 
amounts to 0.1 to 0.6% of the MS durations, depending on the 
mass of the staiQ By that time, the mass fractions of some nu- 
clides have already been substantially modified in their core. 
Most noticeable is 3 He, the evolution of which, at the center 



of the stars, is shown in Fig. 19 (top panel) for selected stellar 



masses. He is produced by deuterium burning during the pre- 
MS phase in all our stars, but is destroyed in more massive, hot- 
ter stars through 3 He + 3 He. The net result, at the time of ZAMS, 
is an enrichment of 3 He, with respect to its initial abundance, in 
the core of stars less massive than 1.1 M Q and a depletion in stars 
more massive than this mass. During the MS phase, the evolu- 
tion of 3 He is dictated by the pp chain. 

The ages of our models are taken relative to the initial mod- 
els. The ages at ZAMS are thus not zero, but reflect the time 
necessary for the initial models to get to thermal equilibrium and 
burn 0.25% of hydrogen in their core (see Sect. 3.2 for further 



discussion). The ZAMS model of our solar track, for example, 
has an age of 190 million years, and it takes an extra 4.36 billion 
years for the ZAMS model to reach solar luminosity, resulting 
in an age of the solar model of 4.56 x 10 9 yr. The TAMS is set at 
the time when X C (H) has decreased to 10~ 5 . 

Three reference points are further identified along the MS 
tracks to locate 'homologous points' that have comparable evo- 



1 This 'pre-MS' time has no physical meaning since the pre-MS 
phase is not followed in our simulations. 
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Fig. 19. Evolution of several quantities in raw models at 
Zinit = 0.04 as a function of the percentage of hydrogen burned 
at the center. From top panel to bottom: logarithm of the 3 He 
mass fraction at the center of the star logX c ( 3 He), logX c ( 12 C), 
core temperature in 10 6 K, logarithm of the stellar radius rel- 
ative to the solar radius and percentage of time elapsed since 
the first raw model relative to the age at TAMS. Models of 0.5 
(continuous lines), 0.7 (long-short dashed lines), 0.9 (long-short- 
short dashed lines), 1.1 (long dashed lines), 1.5 (short dashed 
lines), and 3.5 M G (dotted lines) are displayed, as labeled in the 
top panel. Each marker in the top panel represents a computed 
model. The vertical line indicates the percentage of burned hy- 
drogen taken to define the ZAMS in all our grids. 



lutionary states. Those reference points are important to en- 
sure easy interpolation between homologous points at different 
masses and metallicities. Their definition is arbitrary to some 
extent, and varies from one author to another depending on the 
characteristics of the grids and their expected usage. The central 
abundances of H and He and the morphology of the tracks in the 
HR diagram tend to be used for the identification of those refer- 
ence points (e.g.|Schaller et al.|199 2 ; Bergb usch & Vanden berg 



1992 Ekstrom et al.| 201 l| l. In this work, we adopt a slightly 
modified procedure that takes the morphology of the tracks in 
the log/? - logL diagram into account. This is made necessary 
by the diversity of morphologies of the tracks in the HR dia- 
gram in the mass range considered in our grids. The procedure 
is detailed in Appendix [BJ Two reference points are related to 
the hooks observed in this diagram (as well as in the HR dia- 
gram) for M > 2.5 M Q , and one reference point is set between 
the ZAMS and the first turn off point (see Appendix |B). 




I 000 7 500 7 000 6 500 6 000 5 500 5 000 4 500 4 000 
Teff (K) 

Fig. 20. Basic (black) and interpolated (red) main sequence 
tracks at Zj n i t = 0.014, in the mass range between 0.5 and 2M 
where the morphology of the tracks differs most from one track 
to the next. The interpolated tracks are obtained by interpolating 
in mass the two encompassing basic tracks. The masses of the 
tracks are indicated on the figure except for tracks between 1 . 1 
and 1.3 M Q where tracks are equally spaced with mass steps of 
0.02 M Q . 



Finally, the last model in our tracks is located at the base 
of the RGB, after the star has crossed the HR diagram. Six ref- 
erence points are thus defined: the ZAMS, the three reference 
points in the MS, the TAMS, and the base of the RGB. Then- 
locations in the HR diagram are illustrated in Fig. [T] for the 
Zinit = 0.014 grid. They define five intervals of time, four in the 
MS and one in the post-MS. 

The basic tracks are then constructed by choosing {100,1 00, 
100, 80, 220} models in the five respective time interval^] de- 
fined by the reference points, plus one model at the base of the 
RGB. The models are distributed linearly in age in the first two 
intervals on the MS, in log X C (H) in the last two intervals on the 
MS, and linearly in age in the post-MS interval. Each basic track 
thus contains 601 points. The reason X C (H) cannot be used in all 
MS intervals is that the hydrogen abundance in the core is not 
a monotonically decreasing function of time in all models (see 
Appendix IB). 



4.2. Interpolated tracks 

Evolutionary tracks at {Z; n j t ,M} couples not available in the ba- 
sic tracks are interpolated, on demand, from the basic tracks. 
The interpolation is done in two steps. Tracks are first inter- 
polated in metallicity: the two consecutive metallicities Zj n i t , a 
and Z; n it : b of our grids that encompass the requested Z; n i t (i.e. 
Z; n it,a - Zjnit - Zjnit,b) and the two consecutive masses M. d and 
in the grids that encompass the requested M (i.e. M a < M < M D ) 
are identified and interpolated in Z; n it to construct tracks at (Zi„; t , 
M a ) and (Z; n i t , M D ). These tracks are then interpolated in mass to 
obtain the requested evolutionary track at (Zjnit, M). The inter- 
polation procedure ultimately reduces to interpolating between 
equivalent models from basic tracks, where the interpolation 
quantity is either mass or metallicity. Interpolation of stellar 
quantities between two models is performed in logarithm for the 
metallicity, mass, effective temperature, surface luminosity and 
abundance mass fractions, and linearly for the age. 



2 When a model required in the basic track does not correspond to a 
computed model in the raw track, it is interpolated from the raw data. 
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Fig. 21. Accuracy of the interpolation procedure in mass. Original tracks at a given mass are compared to tracks interpolated between 
the two basic masses encompassing the given mass. The color-coded data represents the difference of the logarithm of the radius 
(left panel), the logarithm of the effective temperature (middle panel), and the logarithm of the age (right panel) of equivalent models 
between the basic and interpolated tracks. All tracks are at = 0.014. 
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interpolated between Z m \ t 



21 



but for the interpolation accuracy in metallicity. Original tracks at Z m i t = 0.014 are compared to tracks 



; 0.01 and 0.02. The color scales of the three panels are kept identical to those in Fig. 21 



To evaluate the accuracy of the interpolation procedure, we 
construct, for each basic track i of mass M, , a track i' interpolated 
at mass M,- from the adjacent tracks i - 1 and i + 1, and compare 
the interpolated track to the basic track. Figure 20 compares the 
basic and interpolated tracks for the MS phase of Zj n i t = 0.014 
models in the mass range where the morphology of the tracks 
varies the most from one track to the next. Three regions are 
seen to be more sensitive to interpolation. The first is for masses 
between 1.6 and 1.7 M Q , where the amplitude of overshooting is 
increased (see Sect. [2]), leading to longer MS tracks. The second 
region concerns masses between 1.18 and 1.26 M , where the in- 
terpolation inaccuracy results from the transition between purely 
radiative MS core models at low masses to fully convective core 
MS models at high masses. This difficulty is partially palliated 
by the higher density of the tracks between 1 . 1 and 1 .3 M in our 
grids. And the third region concerns lower mass stars, especially 
between 0.6 and 0.8 M Q . 



tively, in any of the Zj n j t = 0.014 interpolated models. In practice, 
a track interpolated in our grids has an even better accuracy than 
those results since the interpolation is done between two succes- 
sive tracks of a grid rather than between one every two tracks of 
the grid. 

A similar analysis, but in metallicity, is displayed in Fig. 22 
It shows that our dense grids lead to a very good interpolation 
in metallicity. This good interpolation accuracy holds true at all 
metallicities considered in our grids. 

Finally, we tested how accurately the radius, effective tem- 
perature, and age can be obtained at any given point in the log 7^ 
- log/? diagram without interpolation, but by considering the 
closest basic track in that diagram. The accuracies reach values 
as high as 1 1, 4.3, and 23% for R, T e g, and the age, respectively, 
which are at least three to four times worse than with interpola- 
tion. Using interpolated tracks is thus strongly recommended. 



Fi gure |2 1 1 shows the deviation of the stellar radius, effective ^ <j / S0 _2 rf lines 



temperature and age (all quantities in log) of the models between 
the basic and interpolated tracks. The deviations are seen to be 
less than 1 % in radius and 0.5% in effective temperature for most 
of the interpolated models, and 5% in age, with a maximum devi- 
ation of 2.6, 1 .3, and 7.8 % encountered for R, T e g, and f, respec- 



The change in surface abundances with time in models including 
atomic diffusion (see Sect. 3.4 1 implies that the measured metal- 
licity at the surface of a star may be lower than its initial metal- 
licity. The distinction between surface and initial metallicities in 
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Fig. 23. Iso-Z sur f lines (continuous lines) in the HR diagram at 
0.8, 0.9, and 1.0 M as labeled in the figure. Also shown are 
the evolutionary tracks at Z; n ; t = 0.014 (long-dashed lines) and 
0.020 (short-dashed lines), from which the iso-Z sur f tracks are 
constructed by interpolation between equivalent models. 
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Fig. 24. Continuous line: effective temperature of the models of 
the iso-Z surt lines at 1 M versus the age of the models. For com- 
parison, the evolutionary tracks at Z m i t = 0.014 (long-dashed 
lines) and 0.020 (short-dashed lines) are also shown. The filled 
circles indicate the positions of the reference points on the MS. 



models including diffusion has already been stressed in the liter- 
ature, for example in relation with MS fitting of subdwarfs (e.g. 


Morel & Baglin|1999 Lebreton et al.|1999 1 


Salaris etal.|2000), 


binary stars fitting (e.g. Lebreton et al.|1999 


, and globular clus- 


ter distance and AF/AZ determinations (e.g. 


Salaris et al.|2000). 



We stress in this respect that our solar track at Z m \ t =0.014 cor- 
responds to [Fe/H] init = +0.048 and not [Fe/H] init = 0.00. 

In this case, to determine the age and the mass of a star from 
its observed position in the HR diagram, it may be erroneous to 
use tracks computed with an initial metallicity corresponding to 
that measured at its surface. A way to circumvent this difficulty 
is to use lines in the HR diagram that link models of equal initial 
mass, having different initial metallicities and ages, and showing 
equal surface abundances. We call these lines iso-Z sur f lines. By 
searching the values of the initial mass passing through the ob- 
served positions, it is possible to determine the initial mass, the 
initial metallicity, and the age of the star. 

Figure 23 illustrates iso-Z sm f lines at Z su ,f = 0.014 at three 



stellar masses, as well as the corresponding evolutionary tracks 
at Zi n i t - 0.014 and 0.020 from which they are constructed. The 
ZAMS models of the iso-Z sul f lines are identical to the evolu- 
tionary tracks at Z; n j t = Z sur f . Evolved models on the MS of the 
iso-Z sult lines then deviate from the evolutionary tracks. In the 
post-MS phase, the iso-Z sur f models come back close to the evo- 
lutionary models at Z m i t = Z sur f as the surface chemical elements 
are mixed again with the deep layers due to the deepening of the 
convective envelope. 

Iso-Z sur f lines are, of course, not evolutionary tracks. In par- 
ticular, the ages of the models in an iso-Z sult line do not neces- 
sarily increase with model number. This is illustrated in Fig. 24 
which displays the iso-Z sur f track at 1.0 M and Z sul f = 0.014 
The models in the post-MS phase are seen to have an age that 
decreases towards the end of the track, getting back close to the 
ages of the evolutionary track at Z m n = 0.014. 



4.4. Isoch rones 

Isochrones consist of a succession of models all having the 
same age and initial metallicity, but with increasing initial stellar 




4,05 4,00 3,95 3,90 3,85 3,80 3,75 3,70 3,65 3,60 
log Teff (K) 

Fig. 25. Isochrones in the HR diagram for five different ages as 
labeled in Gyr next to the last model of the isochrones (displayed 
in thick black lines). Also shown in gray are the evolutionary 
tracks for initial masses from 0.8 to 2.6 M as labeled next to the 
ZAMS models of the tracks. The main sequences of the evolu- 
tionary tracks are shown in dark gray and the post-MS phases in 
light gray. Isochrones and tracks are at = 0.014. 



masses from one model to the next in the track. The track dis- 
played by an isochrone in the HR diagram corresponds to the po- 
sition of single-aged stellar populations, such as those obtained 
in single-aged stellar clusters or starburts. 

Examples of Z; n i t = 0.014 isochrones between 0.5 and 8 Gyr 
are shown in Fig. [25" 



They are computed from the basic tracks 
by taking care to adequately reproduce the hooks at the end 
of the MS. The details of an isochrone construction procedure 
are given in Appendix [C] We note that the 4.0 Gyr isochrone 
presents more than two hooks in the HR diagram. This is related 
to the development of the convective core, as explained in the 
appendix. 
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Fig. 26. Masses and radii of detached binary stars from Kraus 
et al. ( 201 1 J > (red points) and Torres et al. (2010) (black points) 
in the mass range below 1.1 M B . Isochrones from our grids are 
shown at 14 (continuous line), 5 (long dashed line) and 1 Gyr 
(short dashed line) forZ; n ; t = 0.014. Isochrones at 14 Gyr are also 
shown for Z^t - 0.006 (upper dotted line) and 0.04 (lower dotted 
line). 



5. Comparison with observations and other grids 

In this section, we compare our models with observations of sev- 



eral binary systems (Sect, 
and with some stellar grids from the literature (Sect. |5.3 i 



3Tb, with two open clusters (Sect.|5.2 



5.1. Comparison with binary systems 

Data from detached binary systems provide unique opportuni- 
ties to test predictions of stellar evolution models. We display 
in Fig.[26]the masses and radii of binary stars from |Torres et al. 
( |2010| l and |Kraus et al.| ( |2~01 1| ) that have masses less than 1.1 M . 
We plot in the same figure isochrones from our grids at several 
ages and metallicities. They encompass well the range of masses 
and radii observed in binary stars for masses above 0.7 M Q . For 
lower mass stars, however, the models predict radii that are big- 
ger than what is measured, a conclusion already stressed by 
|Kraus et al.| ( 201 l| l when comparing their measurements with 
predictions of low-mass star models from Baraffe et al. (1998). 
Kraus et al. ( |201 1[ ) attribute this discrepancy to the possibility 
that short-period systems are "inflated by the influence of the 
close companion, most likely because they are tidally locked into 
very high rotation speeds that enhance activity and inhibit con- 
vection". 



Figures |27 j and |28] display four of the binary systems listed 
2010]) in more details, two with masses be- 

The 



by |Torres etaT 

low 1 .0 M Q (Fig. [27) and two with masses around 1 .4 M Q 



positions of the eight stars are shown in both the mass-radius 
plane (upper panels) and the HR diagram (lower panels). Visual 
isochrone fits to each of the four systems are also shown. 

Figure [27] shows the four stars of binary systems 90 (RW 
Lac +TYC 3629-74 0-1) and 91 (HS Aur + BD+47 1350) of 
Torres et al. ( 2010[ >. According to these authors, the metallic - 
ity of those systems is unknown and their age estimated to be 
around 10 Gyr for both systems. Our isochrones confirm that 
both systems have a similar age, of log(age) =s 10.05. According 
to our models, RW Lac and TYC 3629-740-1 would be born with 
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Fig. 27. Isochrone fits to binaries 90 (R W Lac + TYC 362 9-740- 
1) and 91 (HS Aur + BD+47 1350) of |Torres et aLlpOTO) in the 
mass-radius (upper panel) and HR (lower panel) diagrams. The 
positions of the stars in those diagrams are shown with the error 
bars quoted by those authors. Isochrone fitting is performed vi- 
sually. The numbers next to the isochrones are the initial metal- 
licity Zi n i t and the logarithm of the age in years. 



a subsolar metallicity Z; n j t = 0.010, while the initial metallicity 
of HS Aur and BD+47 1350 was above solar with Z init = 0.018. 
The current value of their surface metallicities, which must have 
decreased during the course of their evolution due to atomic dif- 
fusion, are predicted to be [Fe/H] = -0.16 (Z sur f = 0.0086) and 
[Fe/H] = +0.11 (Z sur f = 0.016) for the stars in each system, re- 
spectively. The difference in [Fe/H] between the two compo- 
nents of each system is about 0.01. 

Figure 



28 



shows the four stars of binary systems 64 
(V570 Per + HD 19457) and 65 (CD Tau + HD 34335) of |Torres| 
( |2010[ ). The authors quote an almost solar metallicity for 
V570 Per ([Fe/H] = +0.01 + 0.03) and over-solar for CD Tau 
([Fe/H] = +0.08 ± 0.15), and log ages of 8.8 and 9.5, respec- 
tively. We get higher metallicity estimates, of [Fe/H] = +0.08 
(Z = 0.015) for V570 Per and [Fe/H] = +0.30 (Z = 0.024) for 
CD Tau. We also get slightly larger ages, with log(age) = 9.09 
and 9.42, respectively. 

A deeper analysis should, however, be done to evaluate the 
range of metallicities and ages predicted by our grids compat- 
ible with the error bars on the masses, radii, 7^, and logL of 
those stars. This is outside the scope of this article and should be 
the object of a separate study. These comparisons are, however, 
quite encouraging and indicate that reasonable solutions can be 
obtained. 



14 



N. Mowlavi et al.: Stellar mass and age determinations 



zu 
1,9 
1,8 
1,7 
1,6 
1,5 
1,4 
1,3 
1,2 



0.024, 9.42/ 


0.015, 9.09/ 


CDJau | T / 
/ 




HD_34335 ~/ / 




/ _/'TV570_Per 




/ J^>ID_19457 









1,20 1,25 1,30 1,35 1,40 1,45 1,50 1,55 1,60 1,65 1,7C 
M (Msun) 




7000 



6750 6 500 6 250 
Teff (K) 



6000 



Fig. 28. Same as Fig. 27 but for binaries 64 (V570 Per + HP 



19457) and 65 (CD Tau + HD 34335) of |Torres et alH ( |20lO) . 
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Fig. 29. Isochrone fit to the open cluster NGC 3532. The data 
points are from Clem et al. ( 201 1 1. The thick line is an isochrone 
at 300 My. 



5.2. Comparison with open clusters 

We also tested our isochrones on two open clusters. The first one 
is NG C 3532. With an age estimated at 300 Myr ( |Clem et al.] 
201 1) , it allows our models to be checked at the higher mass 
range of our grids. The data collected by Clem et al. ( 201 l| l 
are plotted in Fig. [29] for all stars with My below 6 mag. Our 
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M67 



dist. modulus = 9.72 mag, E(B-V) = 0.04 mag . 
Zinit = 0.0157, [Fe/H] = +0.10, age = 4.1 Gyrf 
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Fig. 30. Isochrone fit to the open cluster M67. The data points are 
from Montgomery et al. ( 1993| l. The thick line is an isochrone at 
4.1 Gyr. 




300 My isochrone adopts a distance modulus (m - M)y = 8.54 
(492 pc) and a reddening E(B - V) = 0.04. The conversion from 
{r eff ,logL} to {(B - V),M V } uses the T e g to (B - V) conver- 
sion relations of [Sekiguchi & Fukugita (2000 ) and [Bo ehm- 
|Vitense| (|T981 ) and the bolometric corrections of |Flower| ( |1996] r 
A 310 My isochrone at [Fe/H] = +0.10 nicely fits the MS lower 
boundary of NGC 3532. This slightly over-solar metallicity 
agrees with the [Fe/H] = +0.0+0.1 value adopted by Clem et al. 
d§JTT) . 

The second cluster is M67. Its turn-off mass of about 
1 .28 M makes this cluster particularly attractive for compari- 
son with isochrone predictions. It displays a clear hook in the 
CMD (see Fig. [30]). |VandenBerg et al.| ( |2007| > claims that stel- 
lar models with the new solar composition have difficulties in 
predicting this hook, while [Magic et all] (|20 1 ) stress the sen- 
sitivity of the hook on the constitutive physics of the models. 
In Fig. 



30 



we show a matching 4.1 Gyr isochrone from our 
models. We take (m - M)y =9.72 from Sandquist ( 2004[ > and 
E(B-V) = 0.04 from |Sarajedini et al.| ( (T999) . The adopted metal- 
licity, [Fe/H] = +0.10, agrees with recent abundance determina- 
tions for this cluster o f +0.03+0.07 by |FneI et al.| pOlO) and 
+0.05±0.02(±0.10) by |Pancino et"aT[ l |2010) . 

As for the case of the binaries discussed in Sect. |5.1[ the rela- 
tively good matches obtained here are encouraging, but a deeper 
analysis should be done to explore the range of cluster parame- 
ters that lead to isochrones compatible with their CMDs, a task 
beyond the scope of this article. 



5.3. Comparison with other grids 



We compare in the HR diagram (Fig. 31 1 a selection of our solar 



metallicity tracks with three sets of models from the literature: 



- the older Geneva tracks (Sch aller et al.||T992] ). Those used 
old solar abundances, had an overshooting parameter of 
a ov = 0.20, and did not include diffusion; 
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Fig. 31. Comparison of selected tracks of this work at 
Zinit = 0.014 (solid thick lines) with the Geneva tracks of Schaller 



et al.|(|1992|l at Z = .02 (long-dashed lines), with the tracks of 
Girardi et al. ( 2000[ ) without overshooting at Z = 0.019 (short- 
dashed lines) and with the tracks of Demarque et al.| ( |2004| l at 
Z = 0.02 (solid thin lines). The tracks from the latter three works 
are shown up to the base of the RGB only for clarity. The label 
next to the location of the ZAMS is the stellar mass of the tracks. 



- the tracks by Girardi et al. (2000) without overshooting. 
Those use old solar abundances and did not include diffu- 
sion; 

- the Yonsei-Yale tracks (Demarque et al. 2004). Those use the 
old solar abundances, include overshooting with a m — 0.20 
and use diffusion. It must be noted that their Z = 0.02 models, 



displayed in Fig. 31 have not been calibrated to the Sun. It 
is their models at Z = 0.0 1 8 1 that have been calibrated to the 
Sun, see Yi et al. ( 2001| l. This explains their redder 1.0 M 
track in Fig. 31 compared to the three other 1.0 M Q tracks. 



It is of course difficult to compare in detail the tracks from differ- 
ent authors given the differences in the constitutive physics and 
chemistry. We can nevertheless draw some general conclusions 
by comparison with our grids. 

Above 1.2 M , the main parameter differentiating the tracks 
is the overshooting parameter. The higher value of this parame- 
ter in the old Geneva and in the Yonsei-Yale models with respect 
to our models leads to wider MS tracks. The 1992 Geneva mod- 
els were computed with a ov = 0.20 in order to reproduce the ob- 
served width of the MS in the HR diagram. Since then, rotation 
is known to also contribute to broadening the MS width ( Talon 
et al. 1997). Our nonrotating models with moderate overshooting 
of a v = 0.10, representative of slowly rotating stars, thus leads 
to narrower MS widths. For illustration purposes, we also show 
in Fig. 31 tracks computed with no overshooting by Girardi et al. 
(2000). They show a narrower MS, compatible with the expla- 
nation given here. 

For masses below solar, the redder models displayed by our 
0.8 M Q track compared to the older Geneva models and to the 
ones from Gir ardi et al.| ( |2000| ) mainly come from including of 
atomic diffusion in our models, a process neglected in the two 



other works. The 0.8 M Q track of Demarque et al. ( 2004 1, who do 



include helium diffusion in their calculations, is a 



the tracks of the old Geneva models and of |Girardi et al. (2000 ), 
supporting our results with diffusion. It is however not possible 
to compare our models in detail with those of Demarque et al. 
(2004 ) because of the differences in both the constitutive physics 
and chemical composition. 

6. Data access 

Basic and interpolated tracks, as well as iso-Z sur f lines, will be 
made available for download on our web sit^EI Each downloaded 
file corresponds to either a given metallicity (for basic and inter- 
polated tracks and for isochrones) or a given surface metallicity 
(for iso-Z sllr f lines), and contains as many tracks as the number 
of stellar masses requested by the user. An extract of the file 
is given in Table [2] It contains a file header that indicates the 
number of tracks included in the file, and then lists the differ- 
ent tracks one after the other. For each track, a track header first 
summarizes general information on the track. Those are 

- the type of track (in the example given in Table[2] it is a basic 
track); 

- the metallicity, either Zinit or Z sur f depending on the type of 
track; 

- the age or the stellar mass of the track, depending on whether 
it is an isochrone or not; 

- the number of models in the track; 

- the model indexes of the reference points in the track. The 
reference points define phases in the following way: 

- phase 20: starts at the ZAMS, 

- phase 21: starts at the second (after ZAMS) reference 
point on the MS, 

- phase 22: starts at the third reference point, which corre- 
sponds to the occurrence of the first turn-off point on the 
MS or to the equivalent point if there is no hook, 

- phase 23: starts at the fourth reference point, which cor- 
responds to the occurrence of the second turn-off point 
on the MS, or to the equivalent point if there is no hook, 

- phase 30: starts at the TAMS, 

- phase 31: base of the red giant branch. 

The data of the models are then listed, one line per model. For 
each model, we give, ordered by column: 

1 . the model index; 

2. the phase of the model; 

3. t (yr): the age; 

4. M/M Q : the stellar mass; 

5. log L/L Q : the logarithm of the surface luminosity; 

6. log 7^: the logarithm of the effective temperature; 

7. logp c (g/cm 3 ): the logarithm of the central density; 

8. log T c (K): the logarithm of the central temperature; 

9. Q cc : the mass of the convective core relative to the stellar 
mass; 

10. X C (H): the hydrogen mass fraction at the center; 

11. X c ( 3 He); 

12. X c ( 4 He); 

13. X C ( 12 C); 

14. X C ( 13 C); 

15. X C ( 14 N); 

16. X c ( 16 0); 

17. X c ( n O); 

18. X c ( ls O); 

19. X c ( 20 Ne); 



SO redder than 3 http://obswww.unige.ch/Recherche/evol/-Base-de-donnees- 
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Table 2. Extract of a file content downloadable from our web site. See Sect. Efor more details. 



# Grid of Geneva stellar models 

# Number of tracks : 39 



Basic track 

Zinit 

Minit 

numModels 
Phase 20 
Phase 21 
Phase 22 
Phase 23 
Phase 38 
Phase 31 



0.0140 

0. 5000 

601 
starting at 
starting at 
starting at 
starting at 
starting at 
starting at 



age 
age 
age 
age 
age 
age 



3477464e+09 yr, 
4900246e+10 yr, 
8427257e+10 yr, 
6593657e+10 yr, 
7696842e+10 yr, 
1418004e+ll yr, 



with first model at 

with first model at 100 

with first model at 200 

with first model at 300 

with first model at 380 

with first model at 600 



# phase Age Mass LogL LogTeff 

20 1.3477464e+09 0.500000 -1.357416 3.604120 

1 20 1.7832714e+09 0.500000 -1.356323 3.604018 

2 20 2.2187964e+09 0.500000 -1.355400 3.603940 



LogRhoc LogTc Qcc 

1.861599 6.945235 0.0000 

1.861825 6.944023 0.0000 

1.862834 6.943293 0.0000 



XcH XcHe3 XcHe4 

7.05600e-01 8.13559e-04 2.79558e-01 

7.00524e-01 8.24166e-04 2.84595e-01 

6.95546e-01 8.28330e-04 2.89541e-01 



20. X c ( 22 Ne); 

21. X S (H): the hydrogen mass fraction at the surface; 

22. X s ( 3 He); 

23. X s ( 4 He); 

24. X S ( 12 C); 

25. X S ( 13 C); 

26. X S ( 14 N); 

27. X s ( 16 0); 

28. X s ( n O); 

29. X s ( ls O); 

30. X s ( 20 Ne); 

31. X s ( 22 Ne); 

32. Mi„i t /M : the stellar mass of the initial model of the evolu- 
tionary track from which this model has been computed; 

33. Xiaii. the hydrogen mass fraction of the initial model of the 
evolutionary track from which this model has been com- 
puted; 

34. yj n ; t : same as X m [ t , but for the initial helium mass fraction; 

35. Zi n j t : same as X m i t , but for the initial metallicity; 

36. Z sur f : the surface metallicity; 

37. R/R : the stellar radius; 

38. logp/p Q : the logarithm of the mean stellar density relative to 
the solar mean density, with p = 1.411 g/cm 3 ; 

39. v m a X : the frequency of stellar oscillation at maximum ampli- 
tude, given by Eq. d3J; 

40. A max /A maxQ : the maximum oscillation amplitude relative to 
that of the Sun given by Eq. Q; 

41. Av: the large frequency separation, given by Eq. Q. 



7. Conclusions 

The new grids of stellar models presented in this paper cover 
a dense distribution of masses between 0.5 and 3.5 M Q and 
metallicities Z m i t between 0.006 and 0.04, corresponding to 
[Fe/H] = -0.33 to +0.54. The high density of tracks in mass and 
metallicity allows obtaining reliable interpolated tracks for pairs 
of values (Zj n ; t , M) that have not been computed. We estimated 
from the checks presented in Sect. 4 that the dense grids lead to 
interpolation accuracies lower than one percent for most of the 
models for and T e g, and better than five percent for the age 



of most of the interpolated models (Sect 
required to test them against modern hig 
data. 



4.2[ ). This precision is 
iquality observational 



Including of atomic diffusion in our models with M < 1 . 1 M 
leads to variations in the surface abundances that should be taken 
into account when comparing with observational data of stars 
with measured metallicities. For that purpose, iso-Z sm f lines are 
computed. Due to these variations in the surface abundances, the 
relations linking Z and [Fe/H] (Appendix [A| become inaccurate 
in low-mass stars to about 5 %. 

In the transition mass range 1.15-1.25 M , isochrones may 
display more than two hooks in the HR diagram as a result of 
the development of the convective core in the MS phase. The 
exact morphology of the isochrone may depend on the numerical 
treatment of convection. At solar metallicity, this transition mass 
range corresponds to ages between 3.8 and 4.2 Gyr. 

The basic tracks are available for download at 
http: //obswww.unige .ch/Recherche/evol/-Base-de-donnees- 
An interface will also be added to compute and download 
interpolated tracks, iso-Z sur f lines and isochrones, on demand, 
within the limits of the parameters space of our grids. 

Acknowledgements. We thank the referee, A. Weiss, for his careful reading of 
the paper and his fruitful suggestions. 
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Appendix A: Transformation between Z and [Fe/H] 

A relation between the mass fraction Z of all elements heav- 
ier than helium and the metallicity [Fe/H] = log[X(Fe)/X] - 
log[X(Fe)/X] Q (where the subscript refers to the photospheric 
abundances of the Sun today) can be established if we make 
two assumptions on the chemical evolution of the Galaxy. In 
doing that, we assume to deal with abundances in the interstel- 
lar medium, from which stars are born. We therefore explicitly 
write Z init = 1 - X ink - F init and [Fe/H] init = log[X(Fe)/X] init - 
log[X(Fe)/X] G in order to distinguish initial abundances from 
photospheric abundances of running models, the latter having 
possibly been modified in the course of stellar evolution by 
atomic diffusion. 

The first assumption concerns the galactic evolution of the 
helium mass fraction F; n j t with metallicity Zj n ; t . A linear relation 
is assumed of the form 



Transformation relations between [Fe/H] init and Z; n j t result 



AY 
AZ 



X Z\y 



(A.l) 



where Fo is the primordial helium mass fraction of the uni- 
verse resulting from Big Bang nucleosynthesis and || reflects 
the chemical enrichment of the Galaxy whereby H is globally 
transformed into He and heavier elements. Equation A.l gives 



Xinif = 1 - Fn - 1 + 



AY 
AZ 



Second, we assume a constant X(Fe)/Z ratio throughout the 
evolution of the Galaxy, i.e we neglect enhancement of a ele- 
ments. We then have, for all Z, 



Xinit(Fe) _ X ainit (Fe) 



from assumptions A. 2 and A. 3 
[Fe/H] init 



log 



l-Fo-fl + H)^ 
1-Fp 

1 + £| + ^ X l(H Fe / H ].n. 



log 



(A.4) 



AZ Zq 



With the values adopted in this paper (Fo = 0.248, 
AF/AZ = 1.2857 and Z o /X o =0.0174), Eqs.|A^4|become: 



/0.013 

[Fe/H] imt = -log 0.04 

\ Anit 

= 0.013 
m " 0.04 + 10- [Fe/H] « 



(A.5) 



It is important to keep in mind that Eqs. A.4 and A.5 relate 
to abundances at the time of star formation. In our grids, these 
relations are also valid during the MS and post-MS models for 
stars equal or more massive than 1.1 M Q . For stars less massive 
than 1.1 M Q , however, these relations are no more valid because 
they do not take into account the alteration with time of the sur- 
face abundances due to atomic diffusion. For example, the so- 
lar model in our grids, obtained when the Zi n ; t = 0.014, 1.0 M 
model has reached the solar luminosity and radius at the age 
of 4.57 Gyr, predicts surface abundances equal to X Q = 0.7524, 
F Q = 0.2346, Z = 0.0131, and Z Q /X Q = 0.0174, in agreement 
with observations. With this value for Z n at the surface of the 



Sun today, however, Eq. |A.5 
This is obviously wrong as, 



i predicts a value for [Fe/H] of 0.018. 
yy definition, [Fe/H] = 0.0. The er- 
ror comes from the use of Eq. |A~5] which does not include the ef- 
fect of atomic diffusion during the past 4.57 Gyr evolution of the 
solar model. The effect is seen to amount to ~ 2% on [Fe/H] for 
the Sun, and is expected to amount to at most ~5% for the other 



low-mass models (see Fig. 14 1. Within this precision, Eqs. A.4 
andlA.5lcan be used for all cases. 



Appendix B: Construction of the basic tracks 

B. 1 . Main sequence 

B.1 .1 . Reference points in the HR diagram 

The tracks of several representative masses in the HR dia- 
gram are shown in Fig. B.l They reveal a morphology of the 
MS that depends on the stellar mass. In the low-mass regime 
(M < 0.9 M G ), the morphology is rather simple: both the effec- 
tive temperature and the luminosity are monotonically increas- 
ing with time during the whole MS phase. At higher masses, 
several complexities appear in the MS morphology: 



1. For M > 0.9 M Q , r eff starts to decrease during the MS. The 
point where it starts to decrease is located at the end of the 
MS for the smaller masses, and occurs earlier as we consider 
more massive stars. 7^ decreases right from the ZAMS for 

(A.2) M>1.6M . 

2. Two turn off points appear in the HR diagram for M >1.2 M . 

3. The surface luminosity is a strictly increasing function 
of time during all the MS for M<\AM Q . For masses 
M > 1 .4 M , however, the surface luminosity decreases after 
the second turn off point, reaches a minimum and increases 
again as the star evolves off the MS. 



(A.3) 
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Fig. B.l. Selected Z^t = 0.014 tracks in the HR diagram cover- 
ing the MS and part of the post-MS to illustrate the different 
morphologies. The main sequence is represented in thick lines, 
the filled circle locating its end. 



To better cope with these different morphologies on the MS, 
we shift from the HR diagram to the log/? - logL plane. This 
removes the first complexity because the stellar radius is always 
a strictly increasing function of time from the ZAMS up to the 
first loop, if existent, as illustrated in Fig. B.2 For constructing 
the basic tracks we use the morphologies in the log/? - logL 
plane. 



B.1.2. Reference points in the log/? - logL diagram 



The tracks in the log R - log L plane are shown in Fig. B.2 They 
display the hooks already identified in the HR diagram for stel- 
lar masses above ~1.2 M Q (Sect. [B~TT| . The construction of the 
basic tracks is made in three steps: 

1 . We first identify the second turn off point in the log R - log L 
plane. To this aim, we define the line connecting the loca- 
tions of the Z AMS and TAMS in that plane (dotted line in 
panel a of Fig. B.3 I, and compute the geometrical distance 
of every model on the track to that line. This distance, called 
function 1, is shown in panel b of the figure. The first mini- 
mum of function 1, searched backward from the end of the 
MS, defines the second turning point. It is shown by a filled 
circle in panel b of the figure. 

2. We then identify the first turn off point of the track. We pro- 
ceed in a similar way as above, but work in the log R- func- 
tion 1 plane rather than in the log R - log L plane: we define 
the line connecting the ZAMS to the second turning point 
(dotted line in panel b of Fig. B.3 i, and compute the geomet- 
rical distance of function 1 to that line; this new distance, 
called function 2, is shown in panel c. The location of the 
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Fig. B.2. Same as Fig. B.l but in the log R - log L diagram. 



first turning point is defined as the minimum of that func- 
tion. It is shown by a filled circle in the panel. 

3. Finally, we define a third reference point on the MS track, 
between the ZAMS and the first turning point, based on the 
distance of function 2 to the line connecting the ZAMS and 
the first turning point in the log R- function 2 plane (dotted 
line in panel c of Fig. B.3 I. That new distance is called func- 
tion 3, and the reference point is defined by the minimum of 
that function. It is shown by a filled circle in panel d. 

The three reference points so defined are summarized in panel 



a of Fig. B.3 The procedure is applicable to masses as low as 
1.2 M Q , even if no hook is visible in the log/? - logL plane at 
that mass, as illustrated in Fig. B.4 for the 1.2 M Q case. 

For masses lower than 1.2 M®, the technique fails because 
there is no deflection point on the MS tracks in the log R - log L 
due to the absence of a convective core. For those low mass 
stars, the three reference points are taken at fixed values of X C (H) 
based on the values found at the higher mass stars. Fig. |B.5| 
displays logL as a function of the core hydrogen mass frac- 
tion for stars in the mass range 1.2 to 1.4M , and locates the 
three reference points on the tracks. From these results, the ref- 
erence points are fixed at X C (H) = 0.26, 0.065 and 0.008 for all 
the stars with M < 1.2 M Q at Zj n i t = 0.014. The same values of 
X C (H) are adopted for the tracks at higher metallicities. At lower 
metallicities, different values are adopted: X C (H) = 0.29, 0.070 
and 0.010 at Z init = 0.010, and Z C (H) = 0.33, 0.090 and 0.015 at 
Zinit = 0.006. At Zi n it = 0.006, these fixed values of X C (H) are used 
to define the reference points of the 1.2 M Q track as well. 



B.1.3. Basic MS tracks 

Once the reference points are defined on the MS, basic tracks are 
constructed by distributing a fixed number of models between 
the reference points. The distribution is done regularly in age 
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Fig. B.3. Upper panel: MS track of the Z init = 0.014, 2 M G models 
in the log R-log L plane. Second to lowest panels: functions used 



to identify reference points on the MS track (see Sect. B.1.2 1 
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Fig.B.4. Same as Fig. B.3 but for the 1.2 M @ models. 
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Fig. B.5. log L as a function of decreasing core hydrogen mass 
fraction, in logarithm, for Zj„; t = 0.014 models of, from bottom 
to top, 1.0, 1.1, 1.16, 1.18, 1.2, 1.22, 1.24, 1.26, 1.28, 1.3, 1.4 
and 1.5 M Q . The locations of the three reference points within 
the MS are shown by filled circles on each track. The dashed 
lines denote tracks for which the reference points are defined by 
fixed values of X C (H) (see Sect. |B.1^2] i. 

in the first two intervals (i.e. between ZAMS and the first refer- 
ence point and between the first and second reference points) and 
regularly in logX c (H) between the two remaining intervals (i.e. 
between the second and third reference points and between the 
third reference point and the TAMS). A regular distribution in 
logX c (H) cannot be done before the first turning point because 
the hydrogen mass fraction in the core is not a monotonically 
decreasing function of time for all masses. 

B.2. Post-main sequence 

The construction of the basic tracks during the post-MS phase is 
easier to perform than that of the MS. We locate the base of the 
RGB with a technique similar as the one used to locate the turn 
off points on the MS. For that purpose, we use the angle between 
the X-axis in the log R - log L plane and the line connecting the 
TAMS and the running point on the post-MS track. The base 
of the RGB is located at the point where this angle starts to in- 
crease. We then distribute a fixed number of models between the 
TAMS and that reference point at the base of the RGB, regularly 
spaced in age. 

Appendix C: Isochrone construction 

Isochrones are computed from the basic tracks by taking care 
to reproduce adequately the hooks at the end of the MS. The 
simplest procedure would consist in considering, in turn, each 
evolutionary stage available in the basic tracks (identified by the 
model indexes in those tracks as described in Sect. |4.1) , and in- 
terpolating in the set of models with the same index from the 
different tracks to get a model at the requested age. The result- 
ing isochrone would then consist of stellar models distributed 
according to the evolutionary stages defined in the basic tracks, 
which are, by construction, well sampled around the hooks in 
the HR diagram. 

This procedure is, however, not applicable due to the appear- 
ance of a convective core during the MS at stellar masses around 
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Fig. C.l. Logarithm of the age, in 10 9 yr, of the models in the 
basic tracks as a function of model index for stellar masses be- 
tween 1.0 and 1.4 M as labeled next to the curves. The models 
are taken at Z m n = 0.014. 
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Fig. C.2. Isochrone at 1 Gyr, Z init = 0.014, in the HR diagram, il- 
lustrating the isochrone computation procedure. Basic tracks are 
shown in grey for different stellar masses, continuous lines for 
the MS and dashed lines for the post-MS phase. The isochrone 
is plotted in black, each filled circle representing a stellar model. 
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Fig. C.3. Same as Fig. C.2 but for isochrones at 3.7, 4.0 and 
4.3 Gyr as labeled next to the isochrones. 



second MS reference point, 0.005 M e between the second and 
third reference points, 0.0005 M Q between the third and fourth, 
0.0001 M Q between the fourth reference point and the TAMS, 
and 0.0005 M in the post-MS phase. We then interpolate in each 
track of this dense grid to obtain a model at the requested age, 
and add it to the isochrone. The resulting isochrone obtained at 
1 Gyr is displayed in Fig. |C.2| 

The issue mentioned above related to the development of 
the convective core during the MS does impact the morphol- 
ogy of isochrones that have their MS hooks at stellar masses be- 
tween 1.18 and 1.26 M . AtZj^t = 0.014, it occurs for isochrones 



around 4 Gyr, shown in Fig. C.3 The 4 Gyr isochrone clearly 



displays the signature of the appearance of the convective core in 
the middle of the MS at stellar masses between 1.18 and 1 .20 M Q 
and of the propagation of the convective core to the whole MS 
at stellar masses between 1.24 and 1.26 M Q . These effect of the 
development of the convective core on the morphology of the 
isochrones is limited to ages within 200 million years around 
4.0 Gyr. The isochrones at 3.7 and 4.3 Gyr, for example, are 
seen in Fig. C.3 to display a normal behavior, the former with 
the two hooks characteristic of more massive stars and the latter 
with no hook characteristic of lower mass stars. 



1 .20 M Q . Figure |C.l| plots the age of the models in the Z = 0.014 
basic tracks as a function of model index for masses between 1.0 
and 1.4 M Q . The age at the end of the MS (indexes above 400) is 
seen to monotonically decrease with increasing stellar masses, 
as expected. But this is not the case in the middle of the MS. 
The age of models at index 200, for example, is decreasing with 
increasing mass for all masses except when passing from 1.18 
to 1.20 M and again when passing from 1.24 to 1.26 M Q . The 
former case corresponds to the mass range at which a convective 
core appears in the middle-to-end of the MS, and the latter case 
to the mass range at which the convective core settles during the 
whole MS phase (see Fig.|7]i. 

We therefore proceed differently. At the given metallicity, 
we compute a very dense grid of models by interpolating in the 
basic tracks with a mass step that depends on the evolutionary 
stage. We take mass steps of 0.01 M between the ZAMS and the 
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